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Abstract 

Oscillators  are  affected  by  drifts  (linear  phase  drift,  linear  frequency  drift,  i.e.  quadratic  phase 
drift)  and  different  types  of  noise  according  to  the  power  law  model  of  power  spectral  density 
(from  /-2  to  /+2  frequency  noise,  i.e.  f~A  to  f°  phase  noise).  Generally,  for  long-term  instability 
characterization  (duration  greater  than  one  hour) ,  drift  coefficients  are  estimated  by  using  least 
squares  whereas  noise  levels  are  obtained  from  the  residuals  by  using  variances  (AVAR,  MVAR, 
TVAR,  ...). 

However,  the  low  frequency  noises,  such  as  random  walk  FM,  induce  very  long  term  fluctua¬ 
tions  which  may  be  confused  with  deterministic  drifts.  This  effect,  due  to  the  non-stationarity  of 
these  noises,  depends  on  the  low  cut-off  frequency  which  must  be  introduced  in  order  to  ensure 
power  convergence  for  low  frequencies.  We  calculate  the  standard  deviation  of  ”  artificial”  drifts 
due  to  long-term  random  fluctuations,  versus  the  noise  levels. 

The  first  interest  of  these  results  concerns  the  estimation  of  the  measurement  uncertainty  of 
drift  coefficients  :  knowing  the  noise  levels  of  an  oscillator  we  calculate  the  standard  deviation 
of  the  artificial  drift  coefficient  due  to  these  noises;  thus,  if  a  ’’real”  deterministic  drift  is  iden¬ 
tified  in  the  signal,  its  coefficients  are  determined  plus  or  minus  the  artificial  drift  coefficients. 
The  standard  deviation  of  the  artificial  drift  coefficients  may  be  considered  as  the  measurement 
uncertainty  of  the  deterministic  drift  coefficient. 

The  second  interest  concerns  the  predictability  of  an  oscillator  affected  by  a  deterministic 
drift.  Thus,  the  knowledge  of  the  drift  coefficient  uncertainties  yields  a  criterion  for  quantifying 
the  reliability  of  a  time  error  prediction. 


1  INTRODUCTION 

We  consider  a  sequence  of  frequency  deviation  samples  composed  of  a  deterministic  part,  i.e.  a 
linear  frequency  drift,  and  a  random  part: 

y(tk)  =  C\tk  +  Co  +  e*.  (l) 

An  estimation  by  least  squares  yields  estimates  Co  and  C\  of  the  real  coefficients  Co  and  C\. 
Denoting  the  interpolated  samples  by  y(t),  we  obtain: 

y{tk )  =  C\tk  +  Co-  (2) 

The  residuals  are  defined  as: 

<+  =  y(^fc)  -  y{tk)  (3) 
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1.1  Random  Fluctuations  and  Deterministic  Drifts 

The  instantaneous  frequency  is  defined  from  the  nominal  frequency  and  the  frequency  deviation 
samples  by: 

Vk  =  (1  +  yk)  (4) 

If  the  sequence  yk  is  not  centered,  there  are  two  possibilities: 

•  the  real  nominal  frequency  is  different  from  the  assumed  nominal  frequency:  this  is  a  problem 
of  inaccuracy  of  the  oscillator; 

•  there  are  long-term  random  fluctuations  (with  period  much  longer  than  the  duration  of  the 
sequence)  which  are  seen  as  constant  over  the  sequence[l]. 

The  same  problem  may  occur  with  linear  frequency  drift. 

It  is  impossible  to  distinguish  a  “true”  deterministic  drift  from  a  “false”  random  drift. 


1.2  Statement  of  the  Problem 

The  Power  Spectral  Density  (PSD)  may  be  modelled  as: 

Sy(f)  =  E  h«-fa  (5) 

a=-2 

•  If  no  deterministic  drift  exists,  what  are  the  relationships  between  the  noise  levels  ha  and  the 
estimated  drift  coefficients  Co  and  Ci? 

•  If  a  deterministic  drift  exists,  what  are  the  uncertainties  of  the  estimated  drift  coefficients  C0 
and  Ci? 

•  In  both  cases,  what  is  the  Time  Interval  Error  (TIE)  due  to  an  extrapolation  of  the  linear 
frequency  drift? 


2  LINEAR  REGRESSION 


2.1  Coefficient  Calculation 

We  consider  N  measurements  (t,,  yt):  {(to,  t/o),  *  •  • ,  (tjv- 1,  yjv-i)},  regularly  spaced  with  a  sampling 
period  r0: 

tfc  =  to  +  k.To  (6) 

We  need  to  know  the  coefficient  of  the  linear  model: 

Vk  =  Citfc  +  Cq  +  ek  (7) 

The  most  probable  coefficient  values,  in  the  sense  of  the  least  squares,  are  given  by: 
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2.2  Estimation  of  the  Uncertainties 

From  (8)  and  (9),  it  is  possible  to  calculate  <t2(Cq)  and  <r2(Ci): 


<r2(C0)  = 


4(21V  -  If 


G» 
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(N  -  1  )N2(N  +  l)2r03  C°V  y"  2  *••&) 

(11) 

(52  w)  =  5^X2<&*%> 

«■  j 

(12) 

°2  (I] ti-Vi)  =  To2  52 12  i-MVi-Vi) 

*  j 

(13) 

cov  (52  w,  52  *i  -  w) =  r°  12 12  *•<!«■»*) 

(14) 

where  (  )  denotes  an  average  over  an  infinite  number  of  identical  processes  (ensemble  average) . 


2.3  Correlation  of  the  Samples 
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Table  1:  Correlations  of  the  yk  samples  versus  the  noise  levels  ha.  C  is  the  Euler  constant:  C  & 
0,5772.  Assuming  a  sampling  satisfying  the  Shannon  rule,  the  high  cut-off  frequency  is  fh  =  ^ . 
fi  is  the  low  cut-off  frequency. 


The  PSD  Sy(f)  is  the  Fourier  Transform  of  the  autocorrelation  function.  Thus,  if  no  real  drift  exists 
in  the  sequence,  we  have: 
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v  0 

which  leads  to  the  results  given  in  Table  1. 

2.4  Mean  Value  Subtraction 

Table  1  shows  that,  for  low  frequency  noises  (/-2  and  /-1  FM),  the  correlations  of  the  samples 
depend  on  the  low  cut-off  frequency  /;.  This  cut-off  frequency  must  be  introduced  in  order  to  ensure 
the  power  convergence. 

If  the  inverse  of  the  low  cut-off  frequency  is  much  larger  than  the  duration  of  the  sequence  [f0,  t^], 
the  very  long  term  fluctuations  (period  «  j^)  are  seen  as  a  constant[l]  (see  Figure  1). 

On  the  other  hand,  the  subtraction  of  the  mean  value  of  the  sequence  cancels  the  dependence  on  //. 
Denoting  the  mean  value  of  the  sequence  by  y  and  the  centered  sequence  samples  by  y'k: 

1  N~l 

(i6) 


y'k  =  yk~y  (17) 

The  subtraction  of  the  mean  value  is  equivalent  to  a  correction  of  the  nominal  frequency  by  a  factor 
(1  +  y): 

^«t'o(l  +  y)(i  +  yfe)  (18) 

After  subtraction  of  the  mean  value,  it  follows  that: 
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Thus: 
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Considering  the  new  linear  frequency  drift  model: 

y'k  =  C'ltk  +  C'0  +  e'k 


we  have: 


a  N2  (N  +  l)2r02  ^  U  'y0 


<72(C[)  = 
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v N2(N  -  l)2(iV  +  1)2to4  VZ. 
It  may  be  demonstrated  that  C[  =  C i  and  then  a2( C ()  =  <j2(Ci). 
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2.5  Estimation  of  the  Residuals 

The  differences  between  the  estimated  drift  and  the  y k  samples  are: 

efc  =  Vk  ~  Cltfc  -  Cq 


The  variance  of  the  residuals  is  given  by: 

a2(e)  -  <r2(y)  +  — — — ^Lff2(Ci)  +  tf2(C0) 


JiCov  (E  Sfib.  Co)  ~  JrCov  (J2  *Jb-yjfc,  Ci)  +  (JV  -  1)toCou  (Co,  Cx) 


The  residuals  don’t  depend  on  the  subtraction  of  the  mean  value: 

ek  ~  yk  ~  C,tk  —  Co  =  yk  —  C]tk  —  Cq 


(25) 
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3  RESULTS 
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Table  2:  Standard  deviation  of  the  drift  coefficients  and  of  the  residuals  versus  the  noise  level  ha 
and  the  duration  of  the  sequence  r.  The  high  cut-off  frequency  is  fh  =  and  the  low  cut-off 
frequency  is  f\. 

Thus,  after  measuring  the  ha  noise  levels,  we  may  estimate  the  uncertainties  <j(C0)  and  <r(Ci)  by 
using  Table  2. 

This  table  shows  that  the  subtraction  of  the  mean  value  cancels  the  dependence  of  Co  on  fh-  For 
high  frequency  noises,  a (C0)  remains  very  close  to  <t(Cq).  Moreover,  neither  cr(Ci)  nor  <r(e)  are 
modified  by  this  subtraction. 

3.1  Measurement  Uncertainties  of  Drift  Coefficients 

If  no  real  deterministic  drift  exists,  the  determination  of  the  drift  coefficients  yields: 

-2<r(Co)  <  Co  <  2a(Co)  with  95.5%  confidence 

-2u(Ci)  <  Ci  <  2<r(Ci)  with  95.5%  confidence 
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Thus,  measuring  a  drift  coefficient  C  within  the  interval  [-2<r(C), +2cr(C)]  is  compatible  with  a 
null  drift  hypothesis  (with  a  risk  of  the  second  kind  of  4.5%). 

On  the  other  hand,  if  a  real  deterministic  drift  exists,  the  estimates  C0  and  C\  converge  toward  the 
real  coefficients  Co  and  C\: 

{Co)  =  C0  and  (C1)  =  C1. 

The  uncertainty  domains  of  the  coefficients  Co  and  Ci  are: 

Co  =  Co  ±  2<r(Co)  with  95.5%  confidence 
Ci  =  Ci  ±  2cr(Ci)  with  95.5%  confidence 

3.2  Frequency  and  Time  Error  Prediction 

3.2.1  Frequency  error  prediction 

If  Co  and  C\  are  estimated  over  a  sequence  of  N  samples  (duration  r  =  Ntq),  what  error  results 
from  an  extrapolation  of  the  linear  model  to  tjv  +  T? 

il{tN  +  T)  =  Ci-(tjv  4-  T)  +  Co  (28) 

The  Total  Frequency  Error  (TFE)  may  be  defined  as: 

TFE{T)  =  y(tN  +  T)  -  y{tN  +  T)  (29) 

The  TFE  is  composed  of  a  Deterministic  Frequency  Error  (DFE) : 

DFE(T )  =  (Ci  -  Ci)(tN  +  T)~  (Co  -  Co)  (30) 

plus  a  random  error  (see  Figure  2) : 

TFE(T)  =  DFE(T)  +  yr  {tN  +  T)  (31) 

yr(ti)  is  a  centered  random  variable  without  drift,  with  a  variance  cr2(yr)  =  Ra. 

Thus,  denoting  t'  =  +  T,  we  obtain: 

<: TFE2(T ))  =  <r2(C0)  +  <r2{Ci).t'2  +  cr2(yr) 

-2 Cov  (C0l  t/r(f'))  -  2 Cov  (Ci,  yT (t'))  ■?  +  2 Cov  (C0,  Ca)  .t'  (32) 

Cov  {Co,  yr(t’))  is  the  covariance  between  the  parameter  Co  estimated  over  the  sequence  [to,  tjv]  and 
the  random  sample  yr  at  the  date  t'  =  +  T. 

3.2.2  Time  error  prediction 

If  a  sequence  of  x(t%)  is  known  over  a  duration  r  (from  to  to  f/v  =  to  +  r),  the  Time  Interval  Error 
(TIE)  at  fjv  4-  T  may  be  defined  as  [2,  3]: 

TIE(T)  =  x{tx  +  T)  -  x(tjv)  -TytNiT  (33) 
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with 


1  ftN+T 

t„,T  =  r  y^dt 
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K  u  £  » 


where  y,  is  the  extrapolated  frequency  deviation  at  f;  and  M  is  defined  as  T  =  Mtq. 

N+M— 1 

tie(t)  =  t0  y  (w-w) 

i=N 

Thus,  denoting  M'  =  N  +  M  —  1,  we  obtain: 
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3.2.3  Example  of  f~2  frequency  noise 

In  order  to  use  (32)  the  covariances  Cov(Cq,  yr{t'))  and  Cov{Cy,  yr(t'))  must  be  calculated: 

/^>  \  2(2JV  —  1)  \  6  x— '  x 

(co,yrM‘)  ~  jv(1V  +  1)  _  IV(IV+ 1)  ^ 

_  2(2N-l)*=t  6  ”=* 

N(N  + 1)  “  R%M>  N(N  +  1)  “  •r*m’ 

=  h-2  [-^  ~  x2(tN  +  T) 


For  Cou (Ci,  yr (£')),  we  obtain: 


Cov(Cuyr(t'))  = 


Therefore,  for  an  /  2  frequency  noise,  the  standard  deviation  of  the  TFE  is: 


f  -  [it + Hsr(*  - +  T> 

It  is  interesting  to  notice  that  the  DFE  and  the  Random  Frequency  Error  are  fully  separated: 

J  (TFE2(T))  =  J2<r2(e)  +  <r2(Cx)(tN  +  T) 
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Thus,  if  T  =  0  (interpolation),  the  standard  deviation  of  the  TFE  is  \/2  times  the  standard  deviation 
of  the  residuals,  i.e.  it  is  the  standard  deviation  between  two  residuals. 

Concerning  the  TIE,  from  (36),  (37)  and  (38),  we  obtain: 


sj  (TIE2(T))  = 


15  tjv 


(9  tpj2  +  l3fjvT  +  4  T2) 


(41) 


4  CONCLUSION:  CHOICE  OF  THE  FREQUENCY  MODEL 


What  is  the  physical  meaning  of  the  low  cut-off  frequency  of  an  oscillator?  Is  it  a  real  feature  of 
low  frequency  noises  or  a  mathematical  trick?  In  practice,  it  is  possible  to  avoid  its  use. 

For  an  f~2  frequency  noise,  the  derivative  of  the  frequency  deviation,  the  ageing  z(t),  is  a  white 
noise: 

y(t)  -  f  z(6)d6  (42) 

Jt0 

where  to  is  the  switch-on  date  of  the  oscillator.  In  this  case,  we  have  assumed  that  the  oscillator 
was  syntonized  and  synchronized  at  to-  fi  is  no  longer  necessary,  y(t)  is  a  centered  random  variable 
whose  standard  deviation  increases  with  0: 


1 

0 


(43) 


What  is  the  “real”  frequency  of  the  oscillator  over  T  <C  0:  its  nominal  frequency  or  its  mean 
frequency  over  T? 

The  answer  depends  on  the  frequency  model: 


•  the  use  of  the  power  law  PSD  model  implies  that  the  nominal  frequency  and  the  hQ  noise 
levels  are  time-independent:  they  are  the  constants  of  this  model.  This  model  is  suitable  for 
free-running  oscillators,  e.g.  frequency  standards  involved  in  the  TAI  computation; 

•  the  determination  of  the  nominal  frequency  as  the  mean  frequency  over  a  sequence  of  finite 
duration  implies  that  the  nominal  frequency  is  time-dependent:  this  nominal  frequency  is 
only  valid  over  this  whole  sequence  but  neither  over  a  part  of  this  sequence  nor  over  another 
sequence.  This  model  is  suitable  for  oscillators  used  for  an  experiment  of  well-defined  duration. 
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Figure  2:  Estimation  of  the  drift  over  a  sequence  of  frequency  deviation  altered  by  f~2  FM  noise 
(above).  The  drift  was  estimated  over  the  first  256  samples  (256  sec).  After  this  time,  the  sequence 
moves  away  from  the  estimated  drift.  This  effect  is  more  obvious  in  the  residuals  (below). 


